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GABOR FRAME SETS OF INVARIANCE - A HAMILTONIAN APPROACH 
TO GABOR FRAME DEFORMATIONS 

MARKUS FAULHUBER 


Abstract. In this work we study families of pairs of window functions and lattices which 
lead to Gabor frames which all possess the same frame bounds. To be more precise, for every 
generalized Gaussian g, we will construct an uncountable family of lattices {At} such that each 
pairing of g with some At yields a Gabor frame, and all pairings yield the same frame bounds. 
On the other hand, for each lattice we will find a countable family of generalized Gaussians {gi} 
such that each pairing leaves the frame bounds invariant. Therefore, we are tempted to speak 
about Gabor Frame Sets of Invariance. 


1. Introduction and Notation 

A Gabor frame (or Weyl-Heisenberg frame) for is generated by a (fixed, non-zero) 

window function g £ L^(K‘^) and an index set A C It is denoted by Q{g,A) and consists of 
time-frequency shifted versions of g. 

We denote by A = (x, uj) £ x a point in the time-frequency plane and use the following 
notation for a time-frequency shift by A: 

(1.1) Tr{X)g{t) = MujTx g{t) = e^^^‘^'*g{t — x), x,uj,tGM.‘^. 

The operators involved in Equation are the translation operator 

(1.2) T^g(t) = g(t - x) 
and the modulation operator 

(1.3) M^g{t) = 

The latter one shifts a function in the Fourier or frequency space, hence the name time-frequency 
shift for the composition of the mentioned operators. 

Remark. The point “ • ” denotes the Euclidean inner product of two column vectors, i.e. ui ■ t = 
{u},t) = uj'^t (e.g. in Equation (II.11) 1. Also, we will use the notation x^ = x ■ x, Sx ■ y = x^S^y 
and Sx'^ = x'^S'^x for x, y £ R^ and S a d x d matrix. 

Remark. We note that the translation and modulation operator do not commute in general, in 
fact 

(1.4) 

This formula is closely related to the commutation relations in quantum mechanics. 

Let A be an index set. The time-frequency shifted versions of the window g with respect to A 
are called atoms and the set 

(1.5) g{g, A) = {7r{X)g | A £ A} 

is called a Gabor system. Q{g,A) is called a frame if it fulfills the frame property 

(1-6) A\\f\\l<Y.\^fM^)9)\'<B\\f\\l V/£L2(k<^) 
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for some positive constants 0 < A < < oo called frame constants or frame bounds. The index 

set A C is called a lattice in the time-frequency plane if and only if there exists an invertible 
(non-unique) 2d x 2d matrix S, in the sense that A = The volume of the lattice is defined 

as 


(1.7) 


vol{A) = I det(S')| 


and the density of the lattice is given by 


( 1 . 8 ) 


5(A) = 


1 

vol{A) 


For more details on frames, Gabor frames and time-frequency analysis we refer to 0 , m, m, 

[ID, m- 


Remark. The reader familiar with the topic of time-frequency analysis will note that we restricted 
ourselves to the Hilbert space case. The spaces usually involved when it comes to time-frequency 
analysis are the modulation spaces 

(1.9) MP{R‘‘) = {feS'{R‘^)\\\f\\l,,<oo}, 
where 

( 1 . 10 ) \\frMP= [ [ \{f,M^T,g)\Pdxduj, 

with (non-zero) fixed g € S(R‘^) (see e.g [16], [17]). We note that Mp(R‘^) is independent of the 
choice of g € S(R‘^) and furthermore, M^(R‘^) = L‘^{R‘^). In time-frequency analysis the spaces 
M^(R‘^) and its dual space M°°{R'^) replace in a natural way the Schwartz space ^(K^^) and its 
dual space, the space of tempered distributions 5'(M‘^) usually used in the field of analysis. In time- 
frequency analysis the test or window functions are often assumed to be in Feichtinger’s Algebra 
50(1^^^) = M^(IR‘^), whose dual space Sq{R'^) = M°°{R^) is the natural space of distributions in 
time-frequency analysis. For further reading on the topic of modulation spaces we refer to m, 

[ 16 ], m- 

From Equation (ira we see that the frame bounds A and B depend on the window g and 
the lattice A. Fixing the density of the lattice, a question arising is how far one can deform the 
window or the lattice or both without destroying the frame property. Results in this direction 
are usually called perturbation or deformation results and some are given in i, m. m. [n]- 
Usually, results concerning deformations of Gabor frames do not describe the behavior of the 
frame bounds explicitly, but rather state whether the frame property is kept at all or not. We will 
present some perturbation results where not only the frame property is kept, but also the frame 
bounds. 

Although many of the definitions and well-known, general theorems are stated for L^(R'^), we 
will state our results only for d = I. 

This work is structured as follows. In Section [2] we recall the basic properties of the symplectic 
group, in Section [3] we recall the corresponding properties of the metaplectic group as well as the 
interplay between the two mentioned groups. Finally, in Section |4] we will state and proof the 
main theorem of this work and present some examples. 


2. The Symplectic Group 

As we want to describe a lattice by a matrix one possible way of describing the deformation 
process is by multiplying the generating matrix with another matrix from the left. In particular, 
all our deformations will by carried out by symplectic matrices. Hence, in this section we want 
to recall some basic facts about the symplectic group and its elements, the symplectic matrices. 
They are widely used in Hamiltonian mechanics and also serve as tools in time-frequency analysis 
0, [S], [16]. All results of this section can be found in full generality in de Gosson’s book 0. 
We will also list further references at many points. Although we will only treat the 1-dimensional 
case, most results can be formulated verbatim for d > 1 which is why we formulate them for the 
latter case. 
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2.1. Symplectic Matrices. 

Definition 2.1. A matrix S G GL{2d,'R.) is called symplectic if and only if 

(2.1) SJS'^ = S^JS = J, 

where J = ^ ^ ^, 0 is the d x d zero matrix and I is the d x d identity matrix. J is called 

the standard symplectic matrix. 

Remark. We note that condition m is redundant. Actually, we have 

(2.2) SJS'^ = J^ S^JS = J. 

From (12.111 we conclude that all symplectic matrices S G Sp{2d, R) must have determinant equal 
to ±1. In fact, if 5 G Sp{2d,M.) then det{S) = 1, see [6], [7], [II]. Also, Sp{2d, R) is a subgroup 
of SL{2d,M.) and in the case d = 1 we have Sp{2,M.) = SL{2,M.). In all other cases where d > 1, 
Sp{2d,M.) is a proper subgroup of SL{2d,M.). 

Lemma 2.2. The set of all symplectic matrices forms a group denoted by Sp(2d,M.). 

Proof. Let Si, S 2 G 5'p(2d,R). It follows from Equation (12.11) that the product 5152 G 5p(2d, R). 
Taking the inverse of the double equality in (12.11) and using the fact that J~^ = — J we see that 
5-1 G 5p(2d, R) if 5 G Sp{2d, R). □ 


It is convenient and commonly used write symplectic matrices as block matrices in the following 
form 

(2.3) S=(c D 

where A, B, C, D are d x d matrices. With this notation we have the following formula for the 
inverse of a symplectic matrix 

(2.4) 5-1 = "jf 

In the case d = 1 this reduces to the well-known inversion formula for a matrix 5 belonging to 
5L(2,R), as A, B,C, D G are scalars. 


2.2. Free Symplectic Matrices. We will now introduce the building blocks of the symplectic 
group, the free symplectic matrices and state that any symplectic matrix is the product of these 
building blocks i, m- 


Definition 2.3. We call a symplectic matrix 5 
if B is invertible. 


A B 
C D 


G 5p(2d, R) a. free symplectic matrix 


Definition 2.4. Let P be a symmetric d x d matrix and let L be an invertible d x d matrix. We 
define the following 2d x 2d matrices 


(2.5) 

^=( 

" 0 

(2.6) 

Vp= ( 

^ I 

(2.7) 

Ml= ( 

^ L-1 
. 0 


I 

0 

0 

I 


0 


which all belong to 5p(2d,R). We call them generator matrices for the free symplectic matrices. 


The name generator matrix is justified by the following propositions. 
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Proposition 2.5. With the notation of Definition 
^ ^ can be factored as 


S = 


C D 


we get that any free symplectic matrix 


( 2 . 8 ) 


s' = y_ 


-DB 




A proof is given in or |15| . It makes use of well-known factorization results and properties 
of symplectic matrices. 


2.3. Generating Functions. Following we will point out connections between quadratic forms 
in {x,x') and free symplectic matrices. The motivation comes from Hamiltonian mechanics. We 
want to describe the motion of a particle depending on two variables usually called position {x) 
and momentum (p) which depend on time ft) and are coupled by Hamilton’s equations 


(2.9) 


i{t) = -^H{x{t),p{t)) 
p{t) = --^H{x{f),pft)). 


Here, H{xft),p{t)) is the Hamiltonian or Hamilton function. For more details on Hamiltonian 
mechanics see [2]. 

Given two different positions x and x' of a particle we want to know the initial and final 
momentum p and p' assuming that the motion is linear, meaning we have the linear system 
{x,p) = S(xfp'). This is equivalent to 


( 2 . 10 ) 


X = Ax' -I- Bp' 
p = Cx' + Dp'. 


In order to solve this system of equations for {p,p'), clearly B has to be invertible. 

In the case of time-frequency analysis the proper way to use and interpret Hamiltonian me¬ 
chanics is by replacing position by time and momentum by frequency. The following proposition 
is again formulated in the context of time-frequency analysis. 


Proposition 2.6. Let S 


A B 
C D 


€ Spfld, R) he a free symplectic matrix. Let P,Q be dx d 


symmetric matrices and let L be a d x d invertible matrix, 
(i) Then we have 


( 2 . 11 ) 


(x, oj) = S{x', Lu') 


uj = dxW{x, x'), 
uj' = —dx'W{x,x') 


where W is the quadratic form 

(2.12) W{x, x') = ^DB~^x^ — B~^x ■ x' + ^B~^Ax'^ 

where DB~^ and B~^A are symmetric. 

(a) To every quadratic form 

(2.13) W{x, x') = — Lx ■ x' + 

we can associate the free symplectic matrix 

(2-14) Sw = (^ PL-^Q^L^ PL-^ ) ' 

We call the quadratic form in (12.1311 the generating function of Sw in (12.141) 


Remark. Note the connection between the generating function and the factorization of a free 
symplectic matrix. 

Theorem 2.7. For every S € Sp(2d,M.) there exist two (non-unique) free symplectic matrices Sw 
and Sw' such that S = SwSw' ■ 
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Corollary 2.8. The set of all matrices 

(2.15) {Vp.Ml.J} 

generates the symplectic group Sp{2d, ffi.). 

3. The Metaplectic Group 

The second way to perforin a deformation of a Gabor frames is to perturb the window. We will 
describe this process by letting some unitary operator act on the window. In fact, we will only 
deal with some special operators called metaplectic. 

The metaplectic group and its elements, the metaplectic operators, are widely used in quantum 
mechanics and in time-frequency analysis. There is a close connection to the symplectic group 
and this interplay might be used to solve problems in quantum mechanics once the solution for 
the corresponding classical problem is known [^. In time-frequency analysis this property can 
be used to deform Gabor frames without destroying their frame property and even keeping the 
optimal frame bounds 0, m- 

Again, since there is not much difference between formulating the results for d = 1 and d > 1, 
we state them in full generality, although we will only consider the case d = 1 later on. The results 
can be found in [5] in full detail. 

3.1. The Group Mp(2d,K). 

Definition 3.1. The metaplectic group Mp{2d, M.) is the connected two-fold cover of the symplectic 
group Sp{2d,'M.). Equivalently, we can define Mp(2d, M) by saying that the sequence 

(3.1) 0 ^ Z 2 ^ Mp{2d, R) ^ Sp(2d, R) ^ 0 
is exact. 

Remark. A sequence 

(3.2) Ao —>■ Ai —>■••• —>■ An —>■ An+l 

of morphisms is called exact, if the image of each morphism is equal to the kernel of the next 
morphism 

(3.3) im(Afe_i-)> Afc) = fcer(Afc-)> Afe+i), k = l,...,n. 

We want to use a more constructive approach to define the metaplectic group. 

3.2. Metaplectic Operators and the Quadratic Fourier Transform. The metaplectic group 
is a group of unitary operators on A"(R") 0,0, m. m- Let £ 5(IR.‘^) be a function in the 
Schwartz space. Following de Gosson 0 we define the following operators. 

• The modified Fourier transform J defined by 

(3.4) Jfj(t) = i-‘^^‘^ [ V>(t')e"^"‘'‘'dt'. 

• The linear “chirps” 

(3.5) \Cpi:{t) = e^^^*'^fj{t) 

with P being a real, symmetric d x d matrix. 

• The rescaling operator 

(3.6) = z”VMet(L)|V'(Lt), 

where L is invertible and n is an integer corresponding to a choice of arg{det{L)), to be 
more precise 

(3.7) mr = arg{det{L)) mod 2tt. 

The class modulo 4 of the integer n appearing in the definition of the rescaling operator (13.61) is 
called “Maslov index” 0, 0. 

As in the section on the symplectic group we will associate quadratic forms to metaplectic 
operators and we will also see the interplay between the symplectic and the metaplectic group. 



6 


M. FAULHUBER 


Definition 3.2. Let Sw be the free symplectic matrix 

L-iQ L-i 
pl-^q-l'^ pl-^ 


(3.8) 


Sw = 


associated to the quadratic form W{t,t') = \Pt^ — Lt ■ t' + ^Qt'^ (compare proposition 

equations (12.131) and (12.141) 1. Let the operators J, V-p and Mp^n be defined as in (lOl) . (1331) and 
(13.61) respectively. We call the operator 

(3.9) Ji^ 

the quadratic Fourier transform associated to the free symplectic matrix Sw- 
For ijj G 5(K‘^) we have the explicit formula 

(3.10) ^nfj{t) = i'^-^ V\det{L)\ f V'(i') e"" dt', 

jRd 

where W{t,t') is again the quadratic form as defined in (12.131) and Definition 13.21 

Remark. Although, all statements in this section were formulated for the Schwartz space 5(M‘^), 
they also hold for Feichtinger’s algebra 5o(K.‘^) as well as for the Hilbert space L^(]R'^) [S]. 

Remark. We will frequently drop one or both of the indices W and n and will write S instead of 
Sw and S or Sw instead of Sw,n- When the context allows, we will also use other indices than 
the ones mentioned. 

Remark. As can be seen by formula p.9p a quadratic Fourier transform is a manipulation of a 
(suitable) function by a chirp, a modified Fourier transform, a dilation and another chirp. This is 
the exact same way in which the fractional Fourier transform is described in [T] with an additional 
dilation in between the modified Fourier transform and the second chirp. Hence, the quadratic 
Fourier transform is an extension of the fractional Fourier transform in the sense that the directions 
in the time-frequency plane are scaled by some factor depending on the angle. For more details 
on the fractional Fourier transform see also [9]. 


and the inverse is 


Proposition 3.3. The operators Sw,n extend to unitary operators on 
given by 

(3.11) Sw.n — Sw*.n*-, 

where W*{t, t') = — W{t', t) and n* = d — n. 

The fact that Sw,n is a unitary operator is clear since, Afp „ and J are unitary. Obviously, 
we have 


(3.12) ' = yp, =Mp-q_, 

and the inverse of the modified Fourier transform is given by 

(3.13) 

We note that 


. -1 


/R'i 


ifit') dt' 




L~^ , — n 


= M 


-L^d-- 




(3.14) 
and hence, 

(3.15) 

Definition 3.4. The subgroup of Z/f(L^(M‘^)) generated by the quadratic Fourier transforms Sw,n 
is called the metaplectic group and is denoted by Mp(2d, R). Its elements are called metaplectic 
operators. 


Sw,n = VqJ ^Mp-l_nVp = Sw*,r 
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Remark. To each quadratic form we can actually associate not one but two metaplectic 

operators as, due to (13.71) . Sw,n and Sw,n +2 = —Sw,n are equally good choices. This reflects the 
fact that the metaplectic operators are elements of the two-fold cover of the symplectic group. 

Theorem 3.5. For every S € Mp{2d,M.) there exist two quadratic Fourier transforms Sw,n an-d 
such that S — n.^ ■ 

The factorization in Theorem 13.51 is not unique as the identity operator can always be written 
as S\y^jiS\Y* . 

Corollary 3.6. The set of all operators 
( 3 . 16 ) 

generates the metaplectic group. 


Without further preparation we introduce the natural projection of the metaplectic group 
Mp(2(i, R) onto the symplectic group S'p(2d, R.), which we will denote by For the details we 
refer to [5]. 


Theorem 3.7. The mapping 


(3.17) 


TT^P : Mp{2d, R; 

Sw,r. 


Sp{2d, R) 
Sw 


which to a quadratic Fourier transform associates a free symplectic matrix with generating function 
W, is a surjective group homomorphism. Hence, 


(3.18) 


Mp 


(ss'^ = n^P (s^ TT^P {s'^ . 


and the kernel of tt^p is given by 


(3.19) ker{'K^P) = 

hence, tt^p : Mp{2d, R) i—>■ Sp{2d, R) is a two-fold covering of the symplectic group. 

Definition 3.8. The mapping tt^p in Theorem 13.71 is called the natural projection of Mp(2d, R) 
onto Sp{2d,M.). 


Remark. The natural projections of the metaplectic generator elements are the symplectic gener¬ 
ator elements. 


(3.20) TT^P (Vp) = Vp, TT^P = Ml, tt^p (j) = J. 


4. Gabor Frame Sets of Invariance 

We prepared the machinery of the symplectic and metaplectic group to the extend we need 
it in order to be able to deform Gabor frames without destroying their frame property. We are 
interested in Gabor frame deformations which leave the frame bounds invariant. 

Definition 4.1. Assume Q{g,A) and Q{g',A') are Gabor frames with the same optimal frame 
bounds A and B. We write 

(4.1) g{g,K)'^g{g',K'). 

Theorem 4.2. Let g{g,A) be a Gabor frame with optimal frame bounds A and B. Let S G 
Mp{2d,M.) with projection Tr^P(S) = S G Sp{2d,M.). Then g{Sg,SA) is also a Gabor frame and 
has the same optimal frame bounds A and B. 

A full proof is given in [5] . 

Remark. For any window a phase factor c G C with |c| = 1 is negligible in the sense that g{g,A) 
and g(cg,A) have the same frame bounds as can directly be seen from equation ()1.6p . 
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Theorem 14.21 is a particular case of the notion of Hamiltonian deformation of Gabor frames 
(see [S]). It tells us under which conditions the frame property as well as the optimal frame 
bounds are kept when a Gabor frame suffers some disturbances. This is a very special case, as 
in general neither the optimal frame bounds nor the frame property might be kept under some 
general deformation of the frame. However, there are cases when the frame property might be 
kept without keeping the optimal frame bounds (see m. [18]). This is usually done by either 
deforming the window and fixing the lattice or the other way round. By Theorem l4.2l we know that 
these approaches are equivalent as long as we stick to symplectic and metaplectic deformations. 

What we will see in the following sections is that it is possible to keep both, the frame property 
and the optimal frame bounds under certain lattice deformations, without changing the window. 
This is due to the fact that generalized Hermite functions, including the generalized Gaussians, 
are eigenfunctions with eigenvalues of modulus 1 of certain metaplectic operators. Hence, the 
corresponding symplectic matrix will deform the lattice, while the window can remain unchanged. 


4.1. Lattice Rotations and the Standard Gaussian. From this point on, we will only consider 
the 1-dimensional case. The most popular 1-dimensional window function is probably the standard 
Gaussian gi{t) = . Although, Gabor frames with Gaussian window have been studied 

intensively, we still want to explore and exploit the Gabor family Q {gi , A) with vol (A) < 1 such 
that the frame property (11.61) is fulfilled [50], jn\ . 

One of the simplest manipulations of our Gabor frame is to rotate the lattice and calculate the 
corresponding window. This means that our lattice is deformed by the rotation matrix 


(4.2) 




COST 

— sinr 


sinr 

COST 


and the corresponding deformation of the window is given by the action of the quadratic Fourier 
transform on the window gi. To derive a formula for the resulting window we use Proposition [T5] 
and Equation (I3.10p . 


(4.3) 


St giit) = i 




|sin(r)| 




where n(T) S {0,1, 2, 3} depends on r and the choice of arg ^ysin(r) j and where 

(4.4) Writ, t') = ^ {{t^ + 1'2) cos(t) - 2tt')) . 

2 sm(r) ^ ^ 

This manipulation is meaningful whenever r ^ kir^ k G 1 a. Performing the calculations, we get 

(4.5) S^giit) = = cgi(t), 

with |c| = 1. Hence, we have the result 

(4.6) giguA) = gigi,SrA), 

which means that the frame bounds of a Gabor frame with window gi stay invariant under rotation 
of the lattice. 

Although, the ambiguity function of the standard Gaussian gi = is well-known and, 

though, it is an easy exercise to compute it, we will still do the calculations, as the procedure will 
be used extensively in a somewhat more general form in the rest of this work. For the definition 
and an interpretation of the ambiguity function see Appendix 
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(4.7) 


Ag, ix, uj)=J^g,{t+ |) g, (t - |) dt 

— f 2i/4g—dt 

Jr 

= V2 [ e-2^(‘'+=""/4)g-27ric^i^^ 

Jr 


= V2i 


V2 


e— 


= e-^^ 2-^l'^Tgx [ujlVt^ 

= 6"’^^ 2-1/^ 31 (a;/\/2) 

= g— 

= e-f(-^+-^) 

Here, denotes the Fourier transform which is given by 
(4.8) :Ff{uj)= ! /(t)e-2—‘dt. 

JR'^- 

What we used in the calculations above are a change of variables and the Fourier invariance of 
the standard Gaussian, .F( 3 i) = 31 . For the latter argument see [15], [IB] . 


4.2. Elliptic Deformations and Dilated Gaussians. In Section 14.11 we saw that using the 
standard Gaussian window the Gabor frame bounds stay invariant under a rotation of the lattice. 
We will now generalize this result using ideas from Hamiltonian mechanics. For an introduction 
to Hamiltonian mechanics we refer to [5]. The rotation matrix 


(4.9) 




COST 

— sinr 


sinr 

COST 


determines the flow of the harmonic oscillator with mass m = 1 and resonance ^1 = 1. The 
Hamiltonian of this problem is given by 


(4.10) 


H(x, u), t) 



2 


and Hamilton’s equations are given by 


where A = (x, w)^ and both, x and w depend on r. 

Written in its most general form the Hamiltonian of the harmonic oscillator is given by 

/A TTf \ m.rfix'^ 

(4.12) H{x,uj,t) = — + —-—, 


where m is the mass of the particle and is the resonance. The trajectories will be ellipses in 
standard position with semi-axis ratio mf 2 . 

Assume, we are given the Gabor frame Q{gi,A) with standard Gaussian window and arbitrary 
lattice A, vol{A) < 1 . Any dilation of the lattice by a symplectic matrix can be compensated 
by a metaplectic dilation of the window such that the frame bounds remain unchanged, so 


(4.13) 




We compute that the dilated window is of the form 

(4.14) M^gi(t) = c(2m)^/^e"™‘" = 9mit), 
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where |c| = 1. Next, we compute the ambiguity function Ag^^. 

Agm{x,oj) = V^ [ 

Jr 

(4.15) = ^ dt 

Jr 

Hence, any level set of Agm will be an ellipse in standard position with semi-axis ratio m. It is 
kind of self-evident to examine the behavior of Q[gm,M^^K) under deformations induced by the 
harmonic oscillator given by (14.1211 with fl = 1. The flow is then determined by the symplectic 
matrix 


(4.16) 


Sr.m. — 


COST — sinr 


—msmr cost 

The corresponding metaplectic operator Sr,m. is determined by 
( 4 . 17 ) = 


V Isinrl 


where n(T) G {0,1, 2,3} depends on the parameter r as well as on the choice of arg ^Vsin and 
where 


(4.18) 


Wr,m{t,t') = „ ^ -kt'^)cOST- 2tt')) 

2smT 


is the generating function of ST,m- The verification of the upcoming formulas (14.1911 and (I4.21|l 
can be found in appendix [BJ We have 

(4.19) 9m(i) — 
with |c| = 1. Hence, equation ()4.19p implies that 

(4.20) gigm,A) = Gigm.Sr,mA), 
where Sr,m is defined as in (I4.16L as well as 

(4.21) Agmix,Uj) = A(^Sr,mgmj ix,Uj), 


where Sr,m is defined as given by equation (14.1711 . Summing up the results we get the following 
theorem. 

Theorem 4.3. Let gm(t) = c (2m)^/"‘e“^™*^ with |c| = 1 and let A C K.^ fee a lattice with 
vol{A) < 1 Let 


(4.22) 


St.-TYI - 


cos T — sm T 


—msmr cost 
fee the deformation matrix acting on the lattice, then 


(4.23) 


gmit') — i 




m 


sinr 




9m{t')dt' 


with |c| = 1 hence, 

(4.24) Q{gnY,A)'^g{gm,Sr,mA). 

Furthermore, for the ambiguity function A (^Sr,mgm^ we have 
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Setting m = 1, Theorem 14.31 implies that the Gabor frame bounds of a Gabor frame with 
standard Gaussian window stay invariant under a rotation of the lattice and hence, theorem 14.31 
is a generalization of a result given in [8] . 

Also, in Theorem 14.31 the symplectic geometry of the lattice remains unchanged, whereas the 
Euclidean geometry of the lattice will change in general. In the case where to = 1, meaning that 
we only rotate the lattice, the symplectic as well as the Euclidean geometric properties are kept. 

In order to derive Theorem 14.31 we used a very geometric approach and a clear picture in mind 
about the flow induced by the harmonic oscillator. The crucial ingredient for theorem 14.31 to 
work is that we could easily and explicitly calculate the eigenfunctions of the metaplectic operator 
involved. We note that similar approaches have already been made in [3], characterizing the 
(dilated) Hermite functions as eigenfunctions of certain localization operators. The geometric 
approach to phase space has also been used in m to construct frames consisting of eigenfunctions 
of localization operators. We only stated theorem 14.31 for the dilated Gaussian window, but the 
result holds for all dilated Hermite functions since they are eigenfunctions of the quadratic Fourier 
transform and have eigenvalues of modulus I. 

1 ±i \ 

^1 1 ) ’ that H generates some version of a rotated hexagonal 

lattice of volume 1. Taking the standard Gaussian gi as window function, any rotated version 
of the hexagonal lattice gives the same frame bounds, hence, we may choose one representative 
among all versions of the rotated hexagonal lattice and denote its generator matrix by Hq. Let 
(5 > 1 be the density of the lattice such that G{gi, is a Gabor frame. Then, using theorem 

14.31 we gain a result closely related to the problem of optimal pulse shape design for LOFDM ^M\ . 

^ ^ 75 ' 


Let = 


2 

^/3 


The matrix S = 




integer lattice. If we want 
(4.26) 

for some g G L^(IR.), then 


generates a lattice which is a 45 degrees rotated version of the 




(4.27) 


Ag{x,uj) = e 


So, the quadratic form in the exponent of the ambiguity function describes an ellipse in standard 
position with semi-axis ratio •\/3. Furthermore, we know that there are uncountably many other 
lattice arrangements which together with g lead to the same frame bounds, namely 

(4.28) g - g (^ 5 , -^Sr,mSZ^'^ . 

This result carries over to and extends the results given in [53] . We will also discuss this result in 
Example O 


4.3. Modular Deformations of Gabor Frames. In this section we will deal with discrete 
deformations of Gabor frames. In particular, the objects of interest are taken from the modular 
group which we define as follows. 

Definition 4.4. The modular group Sp{2,Z) consists of all 2-dimensional symplectic matrices 
with integer entries. 


Remark. Usually the modular group F is defined as the group of linear fractional transformations 
on the upper half of the complex plane which have the form 


(4.29) 


az + b 
cz + d' 


with a, 6, c, d G Z and ad — be = 1. For more details on the modular group see |23] . 
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Consider the integer lattice 1? . The action of the modular group leaves 1? invariant, i.e. 
Bl? = T? for B G Sp{2,I). In other words, B is just another choice for a basis of 7?. In 
particular, any B G Sp{2,Z) provides a basis for Taking any symplectic matrix S G S'p(2,R) 
and any basis B G Sp{2,Z) for this implies that 

(4.30) SI? = SBI?. 

We stay with the integer lattice for the beginning. Let 0{g, -^I?\ <5 > 1 be a Gabor frame and 
let 

(4.31) 


B = 


G Sp{2,I). 


The corresponding metaplectic operator is given by 


(4.32) 


where — 





2 -jjc- b ? 0. In general Bg will differ from g by more than 


just a phase factor as we apply a chirp a modified Fourier transform a dilation and again a chirp, 
but the lattice Kps = remains invariant under a modular deformation. Hence, 

( 4 . 33 ) g{g,Ki^s) = Q{Bg,Ki..s)- 

This result can be extended in an obvious way. Let S G Sp{2,R) and let S G Mp{2,R) be the 
corresponding metaplectic operator, then, 

(4.34) g(Sg, SAi^s) = g(SBg, SAj^s)- 

Therefore, given any lattice A = SAjp there are countably many possible windows which lead to 
the same Gabor frame bounds. We sum up the results in the following theorem. 

Theorem 4.5. Let S G Sp{2,'R), B G Sp(2,I) and let S and B be the corresponding metaplectic 
operators. Let Aj^s = with <5 > 1 and let g be a window function. Then 

(4.35) g(Sg, SAi^s) = g(SBg, SAj^s)- 

Remark. Whereas the deformations in the previous section have been derived from a continuous, 
compact group, the deformations in this section are derived from a discrete, non-compact group. 
Continuous deformation groups will in general change the lattice, whereas the window might stay 
invariant under the corresponding deformation. Discrete deformation groups will in general change 
the window, whereas the lattice might stay invariant under the corresponding deformation. 

4.4. Examples for Generalized Gaussians. We will now illustrate our intuitive geometric 
approach by example. We will use different generalized Gaussians and different lattices. 

Example A. We start with an example inspired by [23]. Let 


(4.36) 


Ah = 


cos(7r/6) cos(7r/6) 

^ V sin(7r/6) 


be a version of the hexagonal lattice of density i5 > 1. Since, Ah is radial symmetric, we choose the 
standard Gaussian gi as window function and have the Gabor system g{gi, Ah) which is a Gabor 
frame. We apply the dilation matrix on the lattice and the rescaling operator M 3 - 1/4 on 

the window. By theorem 14.31 we know that 


(4.37) 

Furthermore, we compute 


gigi^An) = g{M^-i/4gi, M^-i/4Ah). 


1 


f cos(7r/4) sin(7r/4) 


(4.38) M^-i/4Ah = 

which is a 45 degrees rotated version of the integer lattice scaled to have density 5 > 1. 


^ Vs \ — sin(7r/4) cos(7r/4) 


1 


Z^ = 

Vs ^ 
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(a) T = ■K(VI 





(c) r = —7r/12 (d) r = —-Kid 


Figure 1. Illustration of the action of ST,m on the lattice and of 5'r,m on the 
ambiguity function. The small ellipses illustrate the ambiguity functions 
centered at lattice points. The ellipses centered at the origin indicate flow lines 

of the harmonic oscillator. 


The ambiguity function of gi is given by 

(4.39) Agi{x, w) = 

and the ambiguity function of M^-i/igi = g^ is given by 

(4.40) 

Applying the matrix 


Agj^{x,u}) = e y, 

✓3 


U3 


(4.41) 
derive 

frame bounds unchanged and we have 


cos T sin r 

_ sin T 

V3 


COST 


derived from the flow of the harmonic oscillator with mass to = -^ on the lattice will leave the 


(4.42) 


1 


gig,,AH)^g{gi,^s^z^]^g[g. 


The deformation process is illustrated in the Figure [TJ 




Vd 


Example B. Let Ai^s = be the scaled integer lattice of density <5 > 1 and let gi{t) = 

2 i/4g-7rF standard Gaussian. The standard symplectic form J belongs to the modular 

group Sp{2,Z). Hence, JAj^s = Aj^s and Jgi = cgi with |c| = 1. In this case neither the change 
of basis, nor the metaplectic operation have an effect on the Gabor frame. This is due to the fact 
that we could interpret the change of basis as a rotation of the time-frequency plane. Using the 
quadratic representation of gi, the ambiguity function, we see that a rotation does not have an 
effect since it reads Agi{x,u}) = \ 
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Let US now consider the case of and A = which is a 90 degrees rotated version 

compared to the window in Example El which can also be interpreted as the deformation of the 
window under a change of basis using J as basis. We know the ambiguity function of Ag^(x, uj) = 
--(V3x^ + — ) 

e ^ V ^/3/ . yve rotate our lattice by the matrix S-^ and apply the corresponding operator 
S- 1 . on the window. Hence, the ambiguity function of the new window becomes 

^S-fg^{x,uj) = Ag^ 

= Ag, (^M-],SZIX) 

^ g-f 

(4.43) ^ g-f (M3_4/4 S,a)^.(M3_4/4S,a) 

A^S|M3^_3/4M3_i/4S.A 

^ g-f A^S_jM3_i/2S^A 


Therefore, the level lines of the ambiguity function will be ellipses rotated by 45 degrees with 
semi-axis ratio equal to and the lattice will be a scaled version of the integer lattice with 
density d > 1. The action of the metaplectic operator can be interpreted in a very natural and 
geometric way as can be seen by the calculations above. 



(a) AS_^g^ (b) AS_^g^ 

a/3 

Figure 2. Contour plots of the ambiguity functions of two possible generalized 
Gaussians which lead to the same frame bounds for the scaled integer lattice of 

density <5 > 1. 

We could also have rotated the lattice A = in the opposite direction by S:il. Then, we 

would have had -^Jl? as our lattice, which leads basically to the same lattice, but with another 
choice of basis. Hence, the window as well as the ambiguity function would have changed. In this 
case, we can interpret the deformation as a rotation, a shearing or a choice of a different basis. 
The interpretation is left to the reader, but the ambiguity function will be 

(4.44) AS-jl ^g^{x,w) = AS-i^g^/^{x,w) = ). 
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The calculations are left to the reader, but the geometric understanding of the lattice deformation 
is sufficient to understand the deformation of the ambiguity function. The level lines describe the 
same ellipses as before only rotated by 90 degrees (see Figure [2]). 

Remark. Instead of the ambiguity function, which is commonly used in time-frequency analysis, 
we could as well have used the Wigner distribution (see appendix!^. The Hamiltonian formulation 
of the harmonic oscillator describes the behaviour of a system in classical mechanics which can 
also be formulated as a problem in quantum mechanics. The Hermite functions are eigenstates 
of the propagator which corresponds to the flow of the classical problem. Hence, the deformation 
of the lattice corresponds to a problem in classical mechanics and the deformation of the window 
corresponds to the same problem formulated in terms of quantum mechanics. Therefore, using 
the Wigner distribution, all results can also be interpreted for problems in classical and quantum 
mechanics. 


A. The Ambiguity Function 

Definition A.l. The ambiguity function of a function / S L^(R'^) is given by 

(A.l) A/(x, u;)= J^j{t+^)f{t- I) e-2—‘ dt. 

In a similar way we define the cross-ambiguity function of two functions f,g € L^(]R‘^) 

(A.2) = l^j(t+^)g{t- I) ‘ dt. 

The cross-ambiguity function and hence, just as well the ambiguity function are closely related 
to the short-time Fourier transform 

(A.3) Vgf{x, uj) = (/, 7r(A)5) = [ f{t)g{t - dt. 

In fact, they only differ by a phase factor and we have 
(A.4) Ag/(a;,cc) = e""--Vg/(x,a;). 

The difference in the phase factor is due to the fact that the translation and modulation operators 
do not commute. 

(A.5) 

therefore, we have 

Vgfix,uj) = {f,M^T,g) 

= {T-x/2M-u]f ,Tx/2g) 

^ g-ix- (M_^/2T_,/2/, M^/2Tx/2g) 

= (7r(-A/2)/,^(A/2)5) 

= e-"“-“Ag/(a:,ca). 

The ambiguity function is somehow a more symmetric time-frequency representation of a signal 
than the short-time Fourier transform. In addition, the ambiguity function of a dilated Gaussian 
is real-valued as we have seen in (|4.15|) . Furthermore, the ambiguity function already determines 
a function up to a phase factor [16) . The usual interpretation of the ambiguity function is that 
it tells how much a function is spread in time and frequency, similar to the interpretation of the 
Wigner distribution in physics. The Wigner distribution is given by 

(A.7) Wf{x, ca) = (^x + f (^x - 0 ‘ dt. 

It is related to the ambiguity function by the symplectic Fourier transform 
(A.8) >V/(A)= J-(A/)(JA), 
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where fit)e dt is the Fourier transform on 


B. Proof of Theorem 14.31 

We will now prove equations (14.231) and (14.251) of theorem 14.31 Since the ambiguity function 
already determines a function up to a factor of modulus 1, equation (14.231) and (14.251) are equivalent. 
First, we note that 

COS T sin r 

\/rn y/m 

—y/msmT y/mcosT 

where Sr = S't,i is a rotation by —r. This means, that imposing the elliptic flow S'r,m on the 
dilated lattice is the same as rotating the lattice by the corresponding angle followed by the same 
dilation. The second ingredient we need in order to perform the proof is the covariance principle 

(B.2) 7r(A)§ = Stt (5-1A) , 

where 7r(A) = is the time-frequency shift as defined in (11.11) . This implies that 

(B.3) A(sf){X)=Af{S-^X). 

This is a classical result m which is also used in [24] and a similar result for the STFT is given 
in [TBj. Hence, using the fact that 

AgiiX) = Ag,ix,u;) = 

A ^Sr^m 9rri^ — A 

= ^5l (('S'r.mM^) A^ 

= e-f (M^Sry^X) 

— g-f A) 

_ g“f 
^ g-f 

= Agm{x,uj). 

This proves equation (14.251) and hence, the proof of theorem 14.31 is complete. 


(B.4) 

we compute 


(B.5) 


(B.l) 


< 5 'r,m — 


= Ad^Sr, 
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